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Qh' The nonlinear dynamics of charged-surface instability development was investigated 

ry^ • for liquid helium far above the critical point. It is found that, if the surface charge 

^^ . completely screens the field above the surface, the equations of three-dimensional (3D) 

potential motion of a fluid are reduced to the well-known equations describing the 3D 

Laplacian growth process. The integrability of these equations in 2D geometry allows 

Ch I the analytic description of the free-surface evolution up to the formation of cuspidal 

,-rt ' singularities at the surface. 

J— H I It is known [1] that the flat electron-charged surface of liquid helium is unstable if the 

^ I electric field strength above (-E+) and inside (-E-) the fluid satisfy inequality 

r-| ■ where g is the free fall acceleration, a is the surface tension coefficient, and p is the fluid 

, ^ ! density. An analysis of the critical behavior of the system suggests that, depending on the 

dimensionless parameter S = {E_'^—E^'^)/Ec^, the nonlinearity leads either to the saturation 
of linear instability or, conversely, to the explosive increase in amplitude. The first situation 
may result in the formation of a stationary perturbed surface relief (hexagons [2] and many- 
electron dimples [3]) in liquid helium. The use of a perturbation theory, with the surface slope 



O . 

Q\ , as a small parameter, allowed the detailed analytic study of such structures in the critical 

^ ! region [4,5]. In the second case, the small-angle approximation fails. The cinematographic 

(^ I study by V.P. Volodin, M.S. Khaikin, and V.S. Edelman [6] has demonstrated that the 

"^ ■ development of surface instability leads to the formation of dimples and their sharpening in 

.^ ■ a finite time. A substantial nonlinearity of this processes calls for a theoretical model that is 

^, free from the requirement for smallness of surface perturbations and adequately describes the 

^ i formation dynamics of a singular surface profile in liquid helium. This work demonstrates 

R-'l that such a model can be developed if the condition E_ ^ E^ is fulffiled, i.e., the field 

above liquid helium is fully screened by the surface electron charge, and if the electric field 

^ ■ far exceeds its critical value, i.e., E^ ^ Ec. 

Let us consider the potential motion of an ideal fluid (liquid helium) , in a region bounded 

by the free surface z = 7]{x,y,t). We assume that the characteristic scale A of surface 

perturbations is much smaller than the fluid depth. We also assume that 

aEZ^ < A < EJ/{gp), 

so that the capillary and gravity effects can be ignored. The electric-field potential ip{x, y, z, t) 
in the medium and the fiuid velocity potential $(x,?/, z,t) satisfy Laplace equations 

vV = o, V2<l> = 0, (1) 
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which should be solved jointly with the conditions at the surface 

87rp$i + 47rp(V$)2 + (Vv?)2 = E_2, z = r]{x,y,t), (2) 

?7t = $^ - V±77 ■ V^'l', z = r]{x,y,t), (3) 

ip = 0, z = r]{x,y,t), (4) 

and conditions at infinity 

ip — * —zE^, z -^ — oo, (5) 

$-^0, z^-oo. (6) 

Let us pass to the dimensionless variables, taking A as a length unit, E_ as a unit of 
electric field strength, and A-Er^(47rp)^'^ as a time unit. It is convenient to rewrite the 
equations of motion of free surface z = ri{x,y,t) in the implicit form (not containing the 7] 
function explicitly). Let us introduce the perturbed potential ip = (p + z decaying at infinity. 
One has at the boundary: ip\z=ri = V- I^ is then straightforward to obtain the following 
relationships: 
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1- V^2 



Y7 ^^^ 



z=ri 



1 — V^2 
z=ri ~ ^ 

which allow one to eliminate the rj function from Eq. (3). The dynamic and kinematic 
boundary conditions (2) and (3) then take the form 

$t - ^. = -(V$) V2 - (V^) V2, z = r]{x, y, t), (7) 

ipt-(^, = -Vif ■ V$, z = ri{x, y, t). (8) 

Let now introduce a pair of auxiliary potentials: 

0(±)(x,i/,z,t) = (^±$)/2. 

With these potentials, the fiuid surface shape can be defined by relationship 

r^=(0(+) + 0(-)) , (9) 

z=ri 

while equations of motion (l)-(6) are reduced to the following symmetric form: 

V20(^) = 0, (10) 

0p) = ±0W T (V0(±))^ z = v{x,y,t), (11) 

0(^) -^0, z^ -cx), (12) 

where boundary conditions (11) are obtained by combining Eqs. (7) and (8) with plus and 
minus sign, respectively. 

It is seen that the equations of motion are split into two systems of equations for the 
potentials 0*^+-' and (l)^~\ which are implicitly related by equation for surface shape (9). 
An essential point is that these equations are compatible with either (p^^^ = or 0*^+^ = 
condition. One can readily see that the first condition corresponds to those solutions 



whose amplitude increases with time, while the second condition corresponds to the decaying 
solutions that are of no interest to us. 

Thus, an analysis of the equations of motion of a charged surface of liquid helium reveals 
a solution increasing with t and corresponding to the 0^"-* = condition or, what is the same, 
to the ip + z = ^ condition (the stability of this branch of solutions is proved below). The 
functional relation between the potentials can be used to eliminate the velocity potential $ 
from initial Eqs. (l)-(6). In the moving system of coordinates {x', y', z'} = {x, y,z — t}, one 
has 

VV = 0, (13) 



V[ = dnV ^l + (V^V)^ z' = v'{x', y', t). (14) 

V^ = 0, z' = r]'{x',y\t) (15) 

ip -^ —z', zl -^ — oo, (16) 

where r]'{x',y',t) = rj — t, and (9„ denotes the normal derivative. These equations explicitly 
describe the motion of a free charged surface z' = ri'{x',y',t). They coincide with the 
equations for the so-called Laplacian growth process, i.e., the phase boundary movement 
with velocity directly proportional to the normal derivative of a certain scalar field {(f in 
our case). Depending on the system, this may be the temperature (Stefan problem in the 
quasi- stationary limit), the electrostatic potential (electrolytic deposition), the pressure (flow 
through a porous medium), etc. 

Note that the boundary movement described by Eqs. (13)-(16) is invariably directed 
inward from the surface. Let 7]' be a single- valued function of variables x' and y' at zero time 
t = 0. Then the inequality r]'{x',y',t) < r]'{x',y',0) holds for t > 0. In the initial notations, 

Vix,y,t) <r]{x,y,0)+t (17) 

for any x and y. This condition can be used to prove the stability of the ascending branch to 
the small perturbations of potential (f)^~\ Clearly, the boundary motion at small 0*^~^ values 
is entirely controlled by the potential 0*^"^^ [one should set 0*-^-' = in Eq. (9)] and, hence, 
obeys Eqs. (13)-(16). The evolution of the 0*-"-' potential is described by Eqs. (10)-(12), 
with the following simple boundary condition in the linear approximation: 

0t"^ = -0l"\ z = r]{x,y,t). 

Let the potential distribution at zero time t = be determined by the expression 

0^"^|i=o = 0o(a:,|/,^), 

where 0o is a harmonic function aX z < ri{x,y, 0) decaying aX z ^ — oo. It is then straight- 
forward to show that the time dynamics of the 0*^"-* potential is given by 

0(-) = (f)Q{x,y,z-t). 

This implies that the singularities of the (p^'' function will drift in the z direction, so that 
they will occur only in the z > ri{x,y,0) + t region. Taking into account inequality (17), 



one finds that the singularities always move away from the boundary z = r]{x,y,t) of liquid 
helium. Consequently, the perturbation 0*^") will relax to zero, as we wished to prove. 

Let us now turn to the analysis of the dynamics of surface instability development in 
liquid helium. In the 2D case (all quantities are taken to be independent of the y variable), 
system of Eqs. (13)-(16) is reduced to the well-known Laplacian growth equation (see, e.g., 
[7] and references therein): 

ImiflU) = 1, if = 0. 

In this expression, f = x' + iz' is a complex function analytical in the lower half-plane of 
the complex variable w = ip — iip and satisfying condition f^wdXw^ip — ioo. Note 
that the ip function is a harmonic conjugate to (/?, while the condition ip = const defines the 
electric field lines in a medium. The Laplacian growth equation is integrable in the sense 
that it allows for the infinite number of partial solutions: 
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f{w) = w - it - i ^ a„ In (w - Wnit)) + i ^ a„ In (tf; - w{t)) 



n=l \n=l 



where a„ are complex constants and Im(w„) > 0. The last term is added in order that 
the condition 77 — > be fulfilled at |a;| -^ 00. One can put Im(w) ^ Im(w„); in this case, 
the influence of this term on the surface evolution can be ignored. The functions w„(t) are 
deflned by the following set of transcendental equations [7] : 
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Wn + it + i^a*^\n {Wn - w^) = Cn, 

m=l 



where C„ are arbitrary complex constants. 

Let us consider the simplest {N = 1) solution to the Laplacian growth equation: 

f{w)=w-it + i\n{w-ir{t)), r{t) -\nr{t) = I + 1^- 1, (18) 

where tc is a real constant and the real function r(t) > 1. The shape of a solitary perturbation 
corresponding to Eqs. (18) is given parametrically by expressions 

z{tp,t) = lny'V^2^r2(t), a;(V',t) = V'-arctan(^/r(t)). (19) 

This solution exists only during a flnite time period and culminates in the formation, at time 
t = tc, of a singularity in the form of a flrst-kind cusp at the fluid surface. Indeed, setting 
r = r{tc) = 1 in Eq. (19), one obtains 

2z = |3a;|2/^ 

in the leading order near the singular point (see also [8]). Note that the electric fleld turns 
to inflnity at the cusp: 

The surface velocity also becomes inflnite in a flnite time: 

Vt = Zt\^=o ~ l/V^c - t. 



It is worth noting that the singular solution in the leading order is also true when the field 
above the surface is screened incompletely. The point is that the requirement that the field 
above the surface be small compared to the field in fiuid is naturally satisfied in the vicinity 
of the singularity. 

Let now discuss the infiuence of the capillary effects. One can readily estimate the surface 
and electrostatic pressures near the surface: 

aR-^ ~ ap^'^EZ^ (t, - t)~\ {d^^f ~ Xp^'^E^ (t, - t)-\ 

Insofar as we assumed that A ^ aEz^-, the capillary forces cannot compete with the elec- 
trostatic ones, so that there is no need to take into account the surface forces at the stage 
of cusp formation. 

In summary, we succeeded in finding a broad class of exact solutions to the equations 
of motion of a charged surface of liquid helium. It is remarkable that the solutions ob- 
tained are not constrained by the condition for the smallness of surface perturbations: the 
model suggested describes the free-surface instability development up to the formation of 
the singularities (cusps) similar to those observed in the experiment [6]. 
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